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ABSTRACT 

Recent theoretical work has pointed out that the transition layer between a 
jet an the medium surrounding it may be more complex than previously thought. 
Under physically realizable conditions, the transverse profile of the Lorentz factor 
in the boundary layer can be non-monotonic, displaying the absolute maximum 
where the flow is faster than at the jet spine, followed by an steep fall off. Like- 
wise, the rest-mass density, reaches an absolute minimum (coincident with the 
maximum in Lorentz factor) and then grows until it reaches the external medium 
value. Such a behaviour is in contrast to the standard monotonic decline of the 
Lorentz factor (from a maximum value at the jet central spine) and the corre- 
sponding increase of the rest-mass density (from the minimum reached at the jet 
core). We study the emission properties of the aforementioned anomalous shear 
layer structures in kiloparsec-scale jets aiming to show observable differences with 
respect to conventional monotonic and smooth boundary layers. 

Subject headings: galaxies: jets: general; X-rays: general — radiation mecha- 
nisms: non-thermal; acceleration of particles — methods: numerical; MHD 



Introduction 



The interaction of extragalactic jets with their environment leads, under rather general 
circumstances, to the stratification of the beam of t he jet in the direction normal to its 



vel ocity. The m orpholo gy of FR I sources (e.g ., M87. lOwen. Hardee fc Cornwelll Il989l : 3C 
31, lLainglll996l : 3C 296, lHardcastle. et al.l 119971 ) has been explained in terms of a jet whose 
dynamics is dominated by a boundary shear layer. In such layers the emissivi ty in radio and 



optical maps peaks and, in some cases (e.g., Mrk 501: ICawthorne. et al.lll993l ). the spectrum 
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is rather flat, suggesting that the acceleration of non-thermal particles takes place right at 
the boundary region. Further evidence in favor of such stratification is provided by radio 



polarization measurements, which indicate tha t, towards the 



is highly ordered and parallel to the jet axis (IPerlman. et al 



et bou ndary, the magnetic field 



1999). 



A radial stratification of l ets in FR II radio galaxies has also been observed (ISwain. Bridle fc Baum 



1998: 


Cawthorne, et al. 


1993 


)• 1 


by 


Swain, Bridle & Baum ( 


1998) 



spine of the jet, suggesting that most of the observed radiation co mes from the jet bound 



ary la yer. On the other hand, the rails of low polarization found by ISwain. Bridle fc Baum 



(119981 ) close the the jet boundary indicate that the magnetic field in the layer is either 
axial or toroidal but not radial. Even at parsec scale, the FR II radio jet of 1055+018 ex- 
hibits a transverse structure consis tent with a spine - shear boundary layer jet morphology 



(lAttridge. Roberts & Wardlelll999l ). 



Also sources which are in the borderline dividing FR I and FR II sources show evidences 
of radial jet stratification. A couple of examples of such sources are 3 C 15 and Cen A. In th e 
first case the jet is generally narrower in the optical than in the radio (IDulwich. et al.ll2007l ). 
and simple spine- sheath model may account for the polariz ation angle difference s seen in the 
optical and radio data. For Cen A, the recent X-ray data ( Kataoka. et alJbood ) can also be 
properly explained by a stratified jet model with a radially decreasing velocity field. 



The aforementioned observational evidence along with some others (e.g.. lChiaberge et al 



20001 ) . suggests that the velocity in some jet boundary layers is smaller (but still relativistic) 
than the velocity of the beam itself. On the other hand, laboratory experiments show that, 
almost unavoidably , turbulent interaction layers may develop as jets propagate into a vis- 
cous medium (e.g., iBrown fc Roshkolll974l ). These laboratory shear layers display a radial 
velocity profile roughly monotonic in which the velocity of the jet core smoothly decreases 
until it vanishes at the external medium. 

From a theoretical point of view, the Reynolds number of the jet flow may reach values 
of ~ 10 23 and, hence, it seems unquestionable that such flows will quickly develop a turbulent 
boundary layer which spreads in to the flow and leads to the entrainment and acceleration 
of ambient gas (IDe Youngl 119931 ). However, the nature and the amount of viscosity in rel- 
ativistic jets is still largely unknown. Most probably, the effective viscosity in the lateral 
interaction between the jet an the external medium is a mixture of the turbule nt eddy vis- 
cosity, magnetic viscosity (based on the finite size of the Larmor radius; see, e.g.. lBaanlll980l ) 
and cosmic-ray viscosity (i.e., the viscos ity originated due to the diff usion of the momentum 
carried by energetic particles; see, e.g.. lEarl. Jokipii &: Morfilll ll988). Whichever the origin 
of the viscosity is, the boundary layer itself results from the nonlinear growth of Kelvin- 
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Helmholtz (KH) instabilities on the jet surface. The fact that relativist ic jets are prone to 



these instabilities was demonstrated theoretically in the pioneer work of iTurland Scheuer 



( 11976( 1 



to KH instabilities i n the vortex-sheet approximation ( e.g. 



19781 ; iHarded Il979l ) . iFerrari. Trussoni fc Zaninettil (119821 ) and 



Ferrari, Trussoni & Zaninetti 


Birkinshawl ( 


1991 


) presented 



the first attempts to study (in the linear regime) the influe nce of a finite t hickness shear layer 



in classical and relativistic jets, respectively. Very recently, IHarded (120071 ) has also considered 
the stability properties (in the linear regime) of spine-sheath magnetized, relativistic jets in 
which the magnetic field is parallel to the flow velocity. The development of KH modes in 
sheared relativistic (slab) flows from the linear (analytic) regime to the non-l inear regime by 



means of relativistic hydrodynamic simulations has been carried forward by iPerucho. et al. 
( 120071 ) . Alternative approaches to the linear analysis of the stability of strati fi ed rela t ivistic 



jets and sheared r elati vistic jets have been p erform ed by lHanasz Soil (119961 ) , lUrpinl (120021 ) 
Aloy. et al.l (120021 ) and iMeliani fc Keppensl (120071 ). respectively 



Further support in favor of radial jet stratification as a consequence of the interaction be- 
tween a jet and its environme nt is provided by t hre e-dimen sional hydrodynamic simulations 
of relativistic, large scale jets ( Aloy et aD 1999a jbl 3^ 2000a ) and two-dimensional simulations 
of relativistic magnetized jets ( Leismann et al. 20051 ). These simulations show that radially 
stratified jets are developed out of initially uniform beams. These numerical works display 
also a rather smooth transition layer with decreasing transverse velocity towards the jet 
boundary and high specific internal energy. It must be noted that the viscosity in such 
simulations is of numerical origin, and mimics only qualitatively actual viscous flows. 

In addition to the existence of environmental reasons to produce a radial flow stratifica- 
tion in astrophysical jets, this effect can also be a natural consequence of the jet formation 
process. Jet launching from accretion disks may directly lead to a certain mass flux profile 
and/or magnetic flux profile within the jet. This possi bility has been verified by means 
of axisymm etric MHD simulations both in the c lassical (IPudridtz. Rogers fc Ouvedl 12006 



Fendtl 120061 ) . as well as in the general relativistic (IMcKinney fc Narayanl 120071 ) regime. 



Aloy &: Rezzollal (120061 ) show that, under physically plausible conditions (see Sect. |2J), 
the lateral structure of the interaction layer of a relativistic flo w can be richer than previously 
thought. This is also the case in relativistic magnetized jets (IMizuno et al.ll2008l ). Due to a 
purely relativistic phenomenon (the conversion of specific enthalpy into bulk Lorentz factor 
across a flow discontinuity), the jet velocity may increase right at the contact discontinuity 
(CD) separating the jet and the external medium. Beyond the CD the velocity decreases 
steeply to zero in the radial direction. The growth of the velocity or, equivalently, of the 
Lorentz factor of the beam is associated to the development of a rarefaction wave that 
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emerges form the CD where the density and the pressure decrease with respect to their 
corresponding values at the jet spine. The profiles exhibited by the hydrodynamic variables 
in these layers are non-monotonic, in contrast with the smooth monotonic profiles that have 
been typically discussed in the literature. Thereby, we will refer to them as anomalous or 
AR shear layers. 

As mentioned above, jet shear layers, represent natural sites for particle acceleration, 
providing high-energy cosmic rays and influencing the dynamics of relat ivistic jets in extra- 
galactic radio sources by forming cosmic-ray cocoons (lOstrowskil 120001 ) . The efficiency of 
the acceleration process in these t urbulent shear layers d epends on the particle mean free 
path and on the velocity structure. IStawarz fc Ostrowskil (2002; SO02 hereafter) performed 
a through study of both the acceleration processes acting at jet boundaries and the resulting 
observational effects. They former authors considered the possibility that turbulent standard 
boundary layers of kiloparsec scale jets may substantially contribute to the radiative output 
of the jet. In the present work we will closely follow SO02 working hypothesis, however we 
will replace the mon otonic shear layer kinem atic structure they assumed by an anomalous 
one as suggested by lAloy &: Rezzollal (120061 ) . The aim being to show whether AR-shear 
layers imprint any distinctive feature in the radiation produced by the jet as compared with 
standard boundary layers. We show in this paper that, indeed, the radiative output of jets 
bound by anomalous shear layers significantly differs from that of ordinary boundary layers. 
Thus observations may confirm or rule out the existence of such anomalous shear layers in 
kiloparsec scale jets and, consistently, this may be used to constrain hardly known physical 
parameters in these regions. 

The plan of the paper is to first show [§ [2] the influence of magnetic fields in the profiles 
of the physical vari a bles i n ano malous boundary lay ers. In this way we extend the results 
of lAloy fc Rezzollal (120061 ) and iMizuno et al.l (120081 ) to the adequate parameter range. In 
Sect. Owe summarize the basic model developed by SO02 and adapt it to account for the 
kinematic differences when applied to anomalous shear layers whose profiles will be discussed 
in Sect. [2j Section H] provides the spectral energy distribution obtained for the different 
models of kinematic shear layers shown in the previous section. We discuss our results in 
Sect. E]and sum them up in Sect. 15.11 



2. Magnetohydrodynamic boosting in relativistic jet boundary layers 

The dynamics of a relativistic magnetohydrodynamic jet in an external medium can 
be treated as the motion of two fluids, one of which (the jet) is much hotter and at higher 
(or equal) pressure than the other one (the external medium), and is moving with a large 
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tangential velocity with respect to the cold, slowly moving fluid. Furthermore, the jet can 
be magnetized while the ex ternal medium magnetic field is much smaller than that of the 
jet. lAloy fc Rezzollal (120061 ) found that if the specific enthalpy of the jet is sufficiently large, 
a net conversion of internal to kinetic energy can be produced through a purely relativistic 
channel. The reason for its exclusively relativistic character is that, the evolution of contact 
or tangential discontinuities (like, e.g., the one that separates laterally the jet from the 
external medium), is governed by a genuine coupling between the specific enthalpy and 
the Lorentz factor at both sides of the discontinuity which does not exist in Newtonian 
(magneto-)hydrodynamics. Such mechanism yields a substantial hydrodynamic boost along 
the boundary layer between the jet and the external medium. Basically, the boost along the 
jet lateral bou ndaries is due to the work done by the external medium on an overpressured, 
relativ istic jet. iMizuno et al.l (120081 ) extended the validity of the results of lAloy fc Rezzolla 
( 120061 ) to the case in which the jet is magnetized and conclude that the presence of a jet 
magnetic field (either poloidal or toroidal) enlarges the boost that might be obtained by 
purely hydrodynamic means. The magnetohydrodynamic boost is larger for jets carrying 
toroidal than poloidal magnetic fields. 

The present analysis is aimed to show what is the qualitative variation o f the Lorentz fac- 
tor an d of the magnetic field across a jet boundary layer. Therefore, following I Aloy &: Rezzolla 
( 120061 ). we model the interaction between the jet and the external medium as a one dimen- 
sional Riemann problem in Cartesian coordinates. Furthermore, we restrict to the case in 
which the magne t ic fiel d is p erpendicular to the fl ow speed in the jet. As demonstrated by 
Aloy fc Rezzollal (120061 ) and IMizuno et al.l (120081 ) this simple model retains the basic fea- 
tures of the phenomenon. Certainly, a more detailed model would require three dimensional 
relativistic magnetohydrodynamic simulations. 

We solve exactly the Riemann problem set by the jump in the conditions between 
the jet and the external medium assuming that the magnet ic field is perpendicu lar to the 
velocity field. For this purpose we use the method devised by lRomero et al.l (120051 ). The two 
uniform "left" (jet) and "right" (external medium) states which set the proposed Riemann 
problems possess different and discontinuous magnetohydrodynamical properties: the rest- 
mass density p, the total pressure p = p mag + p gas (p mag and p gas being the magnetic and 
the gas pressure, respectively), the components of the velocity normal v n and tangential v t 
to the initial discontinuity, and the magnetization (3 := p m ag/p g as- In the following, we will 
use for the quantities in the jet and in the external medium the subscripts "L" and "R" , 
respectively. 

Like in the purely hydrodynamic case, depending on the conditions of the left and right 
states two qualitatively and quantitatively different solutions develop (Fig. [[]). On the one 
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hand, a pattern of waves can be formed, which we indicate as ^SCS^ (Fig. [T] left panel), 
where refers to the shock propagating towards the left sweeping up the jet, £_> to the shock 
moving towards the right crossing the external medium and C to the contact discontinuity 
between the two. On the other hand, the pattern can change for sufficiently large velocities 
parallel to the jet axis (tangential) and in this case the inward-moving shock is replaced by 
a rarefaction wave thus producing a ^JZCS-^ pattern (Fig. [1] right panel). Once again, we 
point out that the change in the wave pattern does only happen in relativistic (magneto- 
hydrodynamics. Only in the relativistic regime the evolution of a Riemann problem depends 
upon the components of the velocity parallel to the initial discontinuity or, equivalently, 
upon the Lorentz factor at both sides of the initial discontinuity. Furthermore, only in the 
relativistic regime there is a coupling between the momentum and energy fluxes that depends 
upon the specific enthalpy of initial left and right states. 

We notice that the variation of the magnetization in the intermediate left state is less 
than a few percent when a JJZCS-^ solution forms (see Fig. [1] right panel). In case a ^SCS-^ 
pattern develops, the magnetization in the shocked left intermediate state /3* grows due to 
the shock compression, but, in any case, (3* £ (2 — 3) x /3 L . 



Although iMizuno et al.l (120081 ) have already explored the influence of a dynamically rele- 
vant magnetic field in the jet on the aforementioned boost, they restricted their study to only 
four cases and to e xtremely hot jets (p L /p L > 10 5 ), which is appropriate for gamma-ray 
burst jets (see, e.g.. lAlov. et al.l l2000bl ; lAloy. Janka fc Miiller 1 120051 ; lAloy fc Obergaulinger 



20071 ; iBirkl et al.l 120071 ) . Here we consider a parameter space more adequate for kiloparsec 
scale jets. More precisely, we consider warm or cold jets with a ratio p ga s,L/p L < 10 2 and a 
magnetization /3 L , which is varied between 0.1 (slightly magnetized jet) and 10 4 (Poynting 
flux dominated jet)0. The left state will have comparable but larger total pressure than the 
external medium p L = 10p R , will be under dense p L /p R = 10~ 5 , with a fixed bulk Lorentz 
factor T L = [1 — [v *) 2 — (t>™) 2 ] -1 / 2 = 10. We choose to fix the external medium (right state) to 
be cold (p R / p R ~ 10~ 4 ), non-magnetized and non-moving. Precisely, we take for the external 
medium p R = 10~ 6 /p cxt c 2 , p R = 10~ 2 p cxt , v™ = t> R = $ L = 0, where p cxt is a normalization 
constant which allows us to be scale-free. We point out that with this parametrization 
we provide physically plausible conditions for both the jet and the external medium. For 
convenience, hereafter we will assume c = p cxt = 1, unless stated otherwise. 

We explore the resulting Riemann solution as we vary the component of the jet velocity 
normal to the shear layeij^] and the jet magnetization /3 L . We point out that models with 



1 In lMizuno et alJ (2008) a maximum value of f3 L — 1.8 is considered in their model MHDB. 
2 For a relativistic jet it is expected that such a component is ~ <C u* ~ 1. 
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(3 > 100 are actually rather cold, since p ga s,L/p L <C 1 (e.g., for /? > 10 4 , p ga s,L/p L — 10 2 ). 
According to the most broadly acce pted view, even MHD-generated jets may bec ome matter 



dominated at kiloparsec scale (e.g.. iBegelman &: Lilll994t iFendt fc Ouyedl 12004 ) . However 



there are others who advocate jet models w hich are basically electromagnetic entities at 
such large scales (e.g., iBlandfordl 12002k 120031 ) . In view of these two possible extremes, we 
have covered the /3-parameter space with jet models which are both matter dominated and 
Poyinting-ffux dominated at kiloparsec scale. In Fig. [2] we show the value of the Lorentz 
factor reached at the state left to the contact discontinuity T*. This intermediate state is 
made out of boosted or deboosted jet matter depending respectively on whether a ^JZCS-^ 
or a <_SCS-> solution develops. 

Figure [2] illustrates how a magnetohydrodynamic boost, where T* > T L , exists below 
a critical value of the normal velocity (t>™ ~ 0.026 in the case considered here) whereby a 
<JZCS-> solution develops. The critical value of is independent of the jet magnetization 
(note the crossing of all the solutions at the point where T* = T L = 10). 

The magnetohydrodynamic boost is larger if the magnetization of the jet is increased: 
while for a poorly magnetized jet (/3 L = 0.1) the boost may increase the bulk Lorentz factor 
of the layer by ~ 10% — 50%, for a strongly magnetized jet (/3 L > 10 2 ) the Lorentz factor 
increase can be > 100%. Increasing the magnetization beyond /3 L > 10 4 does not produce 
larger boosts (the effect saturates). Below (3 L < 0.1 the results are almost indistinguishable 
from the case j3 L = 0.1 and we have decided not to include more lines in the Fig. [2] for the 
sake of readability. 

The bottom line of this parametric study is that for conditions realizable in kiloparsec 
scale jets, one may find an increase of the Lorentz factor in the transition layer between 
the jet and the external medium ranging from (1.5 — 2) x T L , while the magnetization in 
that layer is basically the same as in the jet if a JIZCS^ pattern occurs. Furthermore, since 
only the strength of the ma gnetohydrodynamic boost, but not its existence, depends on 
the topology of the field (see iMizuno et al.ll2008l ). we point out that it also happens in the 
presence of randomly oriented magnetic fields. 



3. The physical model 

Our model follows very closely that devised by SO02. We assume t hat non-thermal 



particles can be accelerated at th e boundary layers of relativistic jets (e.g., I0strowskilll990 
19981 . l2000l ; iRieger fc Duffyll2004j ) . thereby producing high-energy cosmic rays. As mentioned 
in Sect, [fl the main goal of this paper is to replace the prototype shear layer assumed in 
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SO02 by an anomalous shear layer of the type suggested by lAloy &: Rezzollal (120061 ) and 
whose basic features have been outlined in Sect. [21 At the most basic level, this means to 
replace the kinematic shear layer structure assumed by SO02, namely, a jet sheath where 
the bulk Lorentz factor (T) decreases monotonically (Fig. [3^), by a boundary layer where 
the Lorentz factor develops a spike in which it is larger than at the jet core (Fig. [3b). We 
shall assume that the magnetic field is uniform in the whole shear layer, which is compatible 
with the results of our previous section and a reasonable assumption in our simple model. 

The working hypothesis made by SO02 are also valid in our case since, the only difference 
between monotonic and anomalous layers is simply kinematic (see Sect. 13.21) . Thus, following 
SO02, the maximum electron Lorentz factor which can be obtained by the combination of 
acceleration of particles along the shear layer and the radiative cooling of such particles is 

7eq ^4-10 8 V 8 ^ /2 , (1) 

where B^g is the magnetic field in microgauss and Vg = Va/IO^uls -1 is the Alfven speed 
Va in units of 10 8 cms _1 . For kiloparsec scale jets, typical values of these two parameters 
are Vg ~ 1 and B^g ~ 10, which yields 7 eq ~ 10 8 . 

We assume that the resulting spectral energy distribution of the electrons accelerated 
at the layer consists of a low energy power law n e (7) oc 7~ CT , with a = 2, finished with a 
high-energy component modeled as a nearly mono energetic peak at 7 = ^ eq due to the pile- 
up of accelerated particles caused by losses (see Fig. 1 of SO02). In terms of the Dirac delta 
function, 5, at energies above the injection energy, 7 > 70, the considered quasy-stationary, 
averaged spectrum of electrons can be cast in the form 

n e (7) = a T a exp(-7/7 e q) +65(7-7 eq ), (2) 

where a and h are normalization parameters whose values (a ~ 10 _7 cm -3 , b ~ 10 -14 cm -3 ) 
result from assuming equipartition between the magnetic field energy density u B = B 2 We/8 n 
and the energy densities of both electron spectral components: 

f°° 1 
(7mc 2 )a7 -2 exp(-7/7 eq )rf7 = / (7 mc 2 ) b5(j - 7 eq ) cfy = - u B . (3) 

70 -^70 

We point out that in SO02, the power-law part of the electron spectrum was also 
endowed with an exponential cut-off instead of with a Heaviside function in spite of their 
Eqs. 6 and 7. The function was used in their work only to compute analytic estimates, 
not for the num erical integrations necessary to calculate their spectral energy distributions 



(jStawad 120071 ) . 
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Assuming equipartition between the particle and magnetic field energy in relativistic 
shear ed jets has been proved to be a rather solid theoretical assumption (see, e.g., lUrpin 



20061 ). However, in some cases the values of a and b derived from Eq. [5] and based on 



equipa rtition arguments may yield only lower bounds of the true values. iKataoka &: Stawarz 
(120051 ) argue that powerful jets in quasars and FR II objects can be far from the minimum 
energy condition, and the field strength very likely exceeds the equipartition value. On the 
other hand, upper limits to the inverse Compton radiation of the jet in M 87 imposed by 
HESS and HEGRA Cerenkov Teles copes also indicate that the magnetic field cannot be 
weaker than the equipartition value (jStawarz et al.l 120051 ). 



3.1. Radiative processes in boundary layers of kpc-scale jets 



The relevant radiative processes taking place at the boundary layer of a relativistic, 
kpc-scale jet are, on the one hand, the synchrotron (syn) emission of the spectral family 
represented by Eq. [2] and, on the other hand, their inverse Compton (IC) cooling due to the 
interaction with the previously produced synchrotron photons (synchrotron self Compton 
-SSC-) or with seed photons of the cosmic microwave background (external Compton -EC- 
). These radiative processes are the dominant for large scale relativistic (T > 2) jets at 
distances from the galactic nucleus z > 10 kpc. 

Following SO02, we assume that the magnetic field is randomly oriented in the shear 
layeJE Additionally, we neglect synchrotron self-absorption and, therefore, we limit the 
computation of the resulting spectra to frequencies v > 10 10 Hz. 

Restricted to the case of an isotropic electron distribution n e (7) the synchrotron emis- 
sivity averaged over a randomly oriented magnetic field B can be computed from 

V3e 3 B 



isyn ( 



V) 



R 



CiT 



Air 



dry, 



(4) 



where c\ = 3eB / 47rmc and R(x) is a combination of the modified second order Bessel 
functions (ICrusius fc Schlickeiserlll986l ) as 



R(x) 



X K ( X 

Y K4/3 { 2 



K 1/3 



0.3 



x 



K 



4/3 



K 



1/3 



X 



(5) 



The IC photon emissivity in the source fram^l, hi c (e, f2), at a certain photon energy e 



3 Note that anomalous layers can also form if the magnetic field is randomly oriented, see Sect. [21 

4 In our case the source frame is a cylindrical layer of the boundary region where T is assumed to be 



-10- 



(in units of the electron rest mass, m e c 2 , e = hv/m e c 2 ) and scattering direction Q, is given 
by 

hi c (e, Q) — c / dei (b dTLi I c/7 <z> d£l e (1 — (3 cosip) an^ei, f2j) n e (7, fi e ), (6) 



where nj(ej,f2j) is the seed photon number density as a function of the incident photon 
energy ej and the photon direction n e (e e , fl e ) is the electron energy distribution, ip is the 
angle between the electron and the incident photon directions, 7 = (1 — /?) 1//2 is the electron 
Lorentz factor corresponding to a velocity v P = [3c and Q c = (cos -1 fi c , <f) e ) its direction in 



the plasma source frame (e.g., iDermerl Il995l ). The rest-frame emissivity in CGS units is 
j ic (v, ft) = heh ic (e,Q). 

In the calculation of the emissivity due to IC process (SSC or EC), we limit ourselves 
to the Thompson regime, in which the scattering cross section a = a T = 6.65 x 10~ 25 cm 2 is 
independent of the seed photon energy, which is an adequate approximation in the regime 
7Q < 1. 



Following IDermerl (119951 ). the scattered photons will be beamed along the direction of 



the electron motion (i.e. Q ~ Q c ) and will emerge after scattering with an average final 
energy e ~ |7 2 e syn . We assume as SO02 that the synchrotron radiation and the electron 
distribution are isotropic and, therefore, nj(ej, fli) = n syn (e syn ) / 4tt and n e (7, Q e ) = n e (7)/47r. 
Then, in Eq. [6]one can use a = a T 5(Q — Q e )5 (e — (4/3)7 2 e syn ) in order to obtain the following 
expression of h ssc (e, Q): 

c<j f°° f°° ( I \ 

^ssc( e ^) = ^y o de syn J rf7n syn (e syn )n e (7)5 - -7 2 e syn J . (7) 

Using Eqs. El [TJ the relation j ssc (u,Q) = heh ssc (e,Q) and the optically thin synchrotron 
intensity / syn (i / syn) = jsyn^ = (c/i/47r)e syn n syn (e syn ), where / is the emitting region linear size, 
the integral form of the synchrotron self-Compton emissivity J ssc (i / ) is written as 

f°° ( 3u \ 47 2 

jssci") = °T / Jsyn ~j o ^Ml)^, (8) 



V 4 7 2 / 3 



where 7 low = max(7 , |e) is the lower limit imposed by the restriction 7e syn < 1, which 
guaranties that we stay within the Thompson regime. 

The CMB radiation is strongly anisotropic in the source frame which is moving at a 
bulk Lorentz factor T, providing an external source of seed photons which move opposite to 
the jet direction, i.e., tj) ~ 7r. Thereby, (1 — (3cosip) = (1 + fi c ), where Q c = (cos -1 fi e , e ) 



uniform. 
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is the direction of the colliding electron in the comoving frame of the shear layer, similarly 
to what we made for the case of SSC radiation. In the Thomson regime a = a T 5(Q — 
f2 e )<5 [e — 7 2 e cm b(l + /i e )], since the scattered photon energy is e = 7 2 e cm b(l + /i c ), where e cm b 
is the energy of the incident CMB photon in units of the electron rest-mass. Hence, Eq. [H] 
can be cast as 



cct 



47T 



Jo 



OO /"CO 

^Ecmb 



djn cmh (€ cmh )n e (j)5 [e - 7 2 e cmb (l + /J,)] . (9) 
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In order to compute the EC emissivity produced by the interaction of CMB photons with 
the power-law part of the electron distribution, we model the CMB field as a monochromatic 
radiation with photon energ ies equal to the average value of a thermal black body spectrum 
( IBlumenthal &: Gouldl Il970l ) but blue-shifted by a factor T, i.e., e cm b = r(e* mb ) = 2.7 x 
r kT* /m e c 2 , and a number density n cmb = T u* cmh /((e* cmh )m e c 2 ). The value that we adopt 
for the energy density of the CMB radiation in the rest frame of the CMB is u* mh = 4 • 
10~ 13 erg cm -3 . With these approximations, the EC emissivity associated to the power-law 
part of the electron distribution can be computed analytically 

h NV2 r / , \va- 



8tt 



-u 



cmb 



if 7 2 < 



2.7/cT* 
hv 



[(l + /i)r] 3 / 2 exp 



r(l + /i)(e 



)7 2 



cmb/ "eq 



-1/2 



2.7kT*T(l+/i) 



(10) 



For the monoenergetic electron peak component, we use the exact black body spectrum 
transformed to the source frame. In thi s case, the expression of n r ,rp h(£rmh) in the boundary 
layer comoving frame is (see Eq. 2.58 of IBlumenthal fc Gouldl Il970l ) 



^cmb(e C mb) = 87rm e 



2 \ 2 

1Tl e C 6 cm b 

h 



cxp 



Tm e c 2 e cmh 
kT* 



- 1 



(11) 



where k is the Boltzmann constant. Substituting Eq. [TT] and 7^(7) = 6^(7 — 7 eq ) into 
Eqj9l taking into account that j EC (u,Q) = heh EC (e,Q) and the restriction imposed by the 
Thompson regime (e cm b7 eq < 1), one obtains 

-1 



[V] 



2ha T b 



(1 + /x) 2 c 2 7, 6 



if v< 



cxp 



Vhv 



1 + H)le q m e c 2 



1 



V 



h 



The total EC emissivity is the sum of Eqs. [T0land[T2| i.e., 



[V] 



(12) 



(13) 
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3.2. Kinematic model of the shear layer 

We model a large scale jet as a uniform cylindrical spine of radius R c with a Lorentz 
factor T c followed by a transition layer where the Lorentz factor changes linearly with radius 
until r = Rj, where the flow velocity either becomes zero (similar to the SO02 model) or 
reaches a maximum from which it abruptly decreases to zero (Figs. [31 II]). For completeness, 
we also consider a boundary layer where the Lorentz factor is uniform. More explicitly 

T n r<R r 

WJR:-R C 



Ur) = { r w + (r c -r J|4 )#J r c <t<r 3 , (14) 



1 r > Rj 

where we take T c = 10 for our reference models. The subscript i = 1 labels the case of 
a monotonic shear layer (r^j = 1), i = 2 represents the case of an anomalous shear layer 
(rj )2 = 15) and i — 3 denotes the case of a jet where the boundary layer is uniform and 
has the same Lorentz factor as the core (Tj^ = 10). The kinematic idealization of the case 
corresponding to an anomalous layer is a rough prototype of the structure displayed by the 
Lorentz factor when a ^JZCS^ solution develops (Fig. CD). We neglect the variation of the 
magnetic field across the boundary layer because it is rather small (see Sect. H]) and we 
assume that the magnetic field is uniform and equal to 10 fiG. 

In order to properly normalize the emissivity properties of the different jet models, 
we choose that all the shear layers transport the same kinetic power as the monotonically 
decreasing one with a jet radius R- h i = 2kpc and a density and pressure equal to the values 
chosen for the left state of the models shown in Sect. [2j This means that we adjust the 
external jet radius of the other two cases to R h2 = 1.25 kpc and R h? , = 1.38 kpc. Since, 



typically, for extragalactic jets, p cxt ~ 10 2 ' — 10 24 gem 3 (jFerrarilll998l ). the fluxes of mass 
and of magnetic field in the jet axial direction are 



<]> = 3 . i(P ( — ) ( ^L.) 2 G cm" 2 ( 1 :> ) 

1 io- 4 g; vi kpc/ 

D \ 2 / 



For each cylindrical shell in the shear layer of radius r and moving with a bulk Lorentz 
factor T(r), we compute the local source frame emissivities for the radiative processes con- 
sidered in Sect. |3] (Eqs. IH [H] and [T3J . Under the assumption that the radiating electrons are 
distributed uniformly in the boundary layer (for a comoving observer), the emissivity of the 
sub-layer at the distance r from the jet axis can be computed as 

f (u*, 6\ r) = V 2 (r, 9*) j I (v= "\ . , // = f ) , (17) 

Jiy ' ' t \> jjiy £>(r,0*) p l-/3(r)/i*y v ; 



13 



where T>(r, 9*) = [r(r)(l — /3(r) /■**)]" is the Doppler factor associated to a sub-layer moving 
with the Lorentz factor T(r) = [1 — /3 2 (r)] at an angle 9* = cos -1 fi* with respect to the 
line of sight. The subscript / = syn, SSC or EC specifies the radiative process. 

The observed flux density from the boundary layer volume V* appropriately modified 
by the flow-beaming patterns is 

S*(u*,9*) = d- 2 / dV*f{u*,9*,r), (18) 

d being the distance to the observer (we take d = lOMpc unless stated otherwise). dV* = 
27cA z rdr corresponds to a volume element of the shear layer with an observed length A z , 
which we assume to be A z = 1 kpc. The integrals involved in the evaluation of Eq. [18] are 
performed numerically and, in order to compute the SSC contribution to the flux density, 
we estimate the SSC emissivity (Eq. [S]) by taking I = (R^ — R c< i)/ sin^*. 

4. Spectral energy distribution of different models of boundary layers for 

kpc-scale jets 

Our main goal is to outline the observational differences between the emission of large- 
scale jets laterally endowed by different types of shear layer. The observed spectral energy 
distribution of the radiation emitted by different kinematic boundary layers in a prototype 
kpc-scale jet is show in Fig. for two different viewing angles 9* = 1° and 9* = 60°. We 
neglect both the absorption of very high energy (VHE) gamma-rays during the propagation 
to the observer and any contribution from a separate non-thermal electron population ac- 
celerated at shocks happening at the jet core. The choice of a small and large viewing angle 
is motivated by the fact that for small viewing angles, 9* 4 ~ l/T^, the larger Doppler 
boosting provided to the emitted radiation by the larger Lorentz factor in AR layers than in 
monotonic layers makes the total flux density of the former case also slightly larger (Fig. [5] 
left panel). The differences in the total flux density are smaller in the synchrotron dominated 
part of the spectrum than in the EC dominated one. At larger angles (Fig. fright panel), 
the flux density predicted for standard layers is 2 to 5 orders of magnitude larger than for 
uniform or AR layers due to the Doppler hiding of the later two cases. The differences are 
very important in the EC dominated part of the spectrum (i.e., in the gamma-ray regime). 

Both, anomalous and uniform boundary layers generate a rather similar spectral pattern, 
particularly at small viewing angles. As we observe the jet at larger values of 9* (Fig. fright 
panel) the smallest Lorentz factor of the uniform case provides a smaller Doppler hiding 
which results into a larger flux density than in case of AR layers. 
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Considering specifically the flux density produced by anomalous layers (Fig. [6]), we 
note two main differences with the standard case studied by SO02. First, the separation of 
the SSC and EC peaks corresponding to IC emission of the power-law distribution of the 
electrons (Eq. [2]) is larger in the former case. For jets with AR layers the frequency separation 
between the EC1 and SSC1 peaks grows with increasing viewing angle from less than one 
order of magnitude (9* ~ 5°) to more than two orders (9* ~ 90°). In case of jets endowed 
by standard layers (Fig. [7]), the frequency separation between the EC1 and SSC1 peaks 
also grows with the viewing angle, but much less than in the previous case: the frequency 
separation between the two peaks is smaller than a factor 8 at 9* = 90°. If a sufficiently large 
number of TeV photons were detected, this difference could be used to distinguish between 
jets with standard and anomalous layers. However, we repeat here the cautionary note of 
SO02: the spectral energy distributions (SEDs) computed in this paper correspond to the 
most optimistic scenario, with a highly relativistic jet spine and very efficient acceleration 
creating electrons with large Lorentz factors up to 7 eq ~ 10 8 . Either smaller values of r c or 
7 cq would reduce the observed very high energy gamma-ray flux and shift the corresponding 
peaks to lower energies. 

The second most important difference between AR and standard radiating boundary 
layers is that the exponential decay of the synchrotron spectrum is shifted towards smaller 
frequencies (delving deeper into the very soft X-ray regime). This has a big impact in the 
effective X-ray spectral index as we will discuss below. 



As noted by iKomissarovi (119901 ) and confirmed by SO02, extragalactic jets with a 
monotonic boundary layer may show a jet-to-counterjet radio-to-optical brightness ratio 
(S* yQ (9*) / S* yn (ir — 9*)) which is smaller than if the jet were uniform. We find that such an 
assertion has to be modified for jets limited by AR-boundary layers. In the latter case, at 
small viewing angles (namely, 9* £ 5° for the considered parametrization) the jet/counterjet 
asymmetry is ~ 2 — 3 times larger than the one corresponding to a uniform jet with the same 
kinetic power (Fig. [8]) . At larger viewing angles, the brightness asymmetry turns out to be 
almost indistinguishable from the case of an equivalent uniform jet (see in Fig. [8] how the 
solid and dashed lines practically overlap each other for 9* £ 20° if T c = 10 or for 9* £ 50° in 
case T c = 3). Certainly, since the brightness asymmetry is a result of the Doppler boosting 
(hiding) of the jet (counterjet), the smaller is the value of the maximum Lorentz factor in 
the anomalous shear layer (F^), the smaller will be the maximum asymmetry at small 9*. If 
the Lorentz factor of the jet spine is small, the jet-to-counterjet brightness ratio reduces in 
absolute value (see the grey lines corresponding to T c = 3, in Fig. M)u However, even for the 



5 The relative increase of the Lorentz factor in an anomalous shear layer is rather independent on the 
Lorentz factor of the jet (left state). Therefore, since we have chosen a parametrization for our prototype 
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relatively small value of T c = 3 the jet-to-counterjet brightness ratio of the anomalous shear 
layer case is larger than the one corresponding to the standard case at T c = 10 independent 
of the viewing angle (e.g., at 9* = 0°, log [S* (OVS* (7r)] ~ 4.35 for the AR model, while 
log [*5syn(0)/S'syn( 7r )] ~ 4.25). If there were a handle on the value of the Lorentz factor 
of the jet, this difference might allow us to discriminate observationally between jets with 
anomalous and standard shear layers. However, given the similarity of the jet-to-counterjet 
brightness ratio between jets with anomalous layers and jets with top-hat profiles (uniform), 
it might be practically impossible to disentangle from this unique value whether the jet is 
shielded by an anomalous layer or by no layer at all. 

Considering the results obtained for standard boundary layers, SO02 concluded that 
the mildly relativistic velocities inferred from the observed brightness asymmetries at tens 
of kiloparsec scales, may correspond to a slower boundary layer and not necessarily to the 
faster jet core, which might be highly relativistic at the observed distances. If anomalous 
shear layers may happen in actual jets, slower spines are preferred and relatively moderate 
maximum shear layer Lorentz factors might also be invoked to explain the observed jet-to- 
counterjet asymmetries. 

In our models, the X-ray radiated power of kpc-scale jets bounded by standard layers is 
dominated by the high-energy electron bump spectral component, implying that the resulting 
synchrotron X-ray flux lies above the extrapolated radio-to-optical continuum (see SO02 and 
Figs. El [7]). This is also true for jets endowed by anomalous boundary layers (see Fig. [H]at 
v* ~ 10 17 Hz). For both types of layers, the spectral slope at X-ray frequencies may be 
rather different from that of the power-law at radio-to-optical frequencies, and the difference 
depends also on the viewing angle. This implies that the effective X-ray spectral index 
«x,efr(#*) computed between hv* = 1 keV and hv% = 5keV as 

can be significantly distinct. We point out that S^(u*,9*) is computed in this paper as the 
value of the flux density calculated using Eq. [181 instead of SO02 who only take S^(u*, 9*) oc 
/ rdrV 2 (r, 6*)R(xl q ), with -R(x* q ) given by Eq. El and x* cq = u*/ci^ q V(r, 9*). 

For highly relativistic jet cores (r c = 10), Fig. [9] (black lines) shows that the most 
noticeable difference between jets with anomalous or uniform shear layers and jets with 
standard limiting boundaries, is the fact that the effective X-ray spectral index of the jet 
a Xefr i s l ar g er than the corresponding to the counterjet a^efi f° r inclinations 9* £ 65°. 



jet in which T$ t 2/T c = 1.5, for the case of T c = 3 we take Tj.2 = 4.5. 
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For standard boundary layer jets ft J Xeff < o^xeSi ^* hcWs (i.e., the counter jet has a steeper 
X-ray continuum as compared to the jet spectrum). Furthermore, jets with uniform and 
AR layers show a much largei@ jet and counter jet effective X-ray spectral index than jets 



with standard layers if the inclination is 8* £ 30°. Indeed 



a- 



j 

X.cfT 



in the former two types of 



boundary layers is even larger than the ft x eff corresponding to jets with monotonic layers 
for viewing angles 8* ^ 30°. 



The large values of the effective X-ray spectral index in jets with uniform and AR layers 
is due to the fact that the decay of the spectrum after the synchrotron peak crosses the X-ray 
band and moves towards optical frequencies for inclinations 8* ^ 60° (Fig. [6]). Consistently, 
the flux density at 5keV decays abruptly and, thus ftx,eff grows. 

In case of jets with more moderate spine velocities (r c = 3; Fig. [9j grey lines), there is 
no viewing angle for which ft J Xeff > «xeff> independent on the shear layer model. However, 
it remains true that the counterjet X-ray spectral index of jets bearing uniform or AR layers 
is much larger than the one corresponding to standard boundary layer jets (there is a factor 
of 2 to 4 difference between them). 

Taking together the results for models with jet spines flowing at T c = 3 and 10, it 
turns out that the range of variability of both a 3 x eS and ft x efT is larger for jets flanked by 
uniform and AR layers than for jets with standard boundary layers. Interestingly, the large- 
scale X-ray emission observed from several radio-loud AGNs exhibits rather different spectral 
characteristics. For inst ance, X-ray spectra of the known jets in quasa rs are very flat (e.g., 



ax ~ 0.23 for 3C 207, Brunetti et al.ll200ll ; a x ~ 0.5 for PKS 1127, ISiemirinowska et al. 



2002 



ftx ~ 0.8 for 3C 273 and PKS 0637). The low effective spectral index is consistent 
with small jet inclinations (c.f. Fig. [HI inset) independent of the shear layer model adopted. 
On the other ha nd, jets in radio galax ies tend to d isplay relatively ste ep X-ray spectra, 
ftx > 1.0 - 1.5 JHardcastle et all l200ll for 3C 66B, IWorrall et all I2OO1I for B2 0206 and 
B2 0755). These values correspond to large inclination angles on Fig. [9J and/or lower jet 
Lorentz factors for jets with standard boundary layers. But, if uniform or AR emitting 
layers are considered, the viewing angles needed to account for such values of the X-ray 
spectral index are smaller (8* e [20°, 25°] if T c = 10, or 8* 6 [40°, 55°] if T c = 3). Very 



recently, lHardcastle et all (120071 ) have reported values of the X-ray spectral index ftx > 2.5 
at distances of ~ 4kpc from the nucleus in Cen A. These large values of ftx are not within 
reach of jet models flanked by standard layers (the maximum of ft Xeff is 1.84 at 8* = 90° 
for the model with r c = 10). Nonetheless, jets with uniform or AR boundaries may yield 
ft J Xeff > 2.5 for 8* > 32° (8* > 59°) if T c = 10 (r c = 3). Indeed, considering that the 



Particularly, in the range 9* £ [60°, 90°], a J x cff (Uniform ar) is 6 to 7 times larger than a 



X, off, standard* 
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inclination of Cen A is ~ 60° — 70° (jSkibo. Dermer fc Kinzerlll994j ). our results suggest that, 
if a uniform or an AR emitting boundary layer were responsible for the emission of Cen A 
at distances from the nucleus £ 4kpc, small values of the spine Lorentz factor (r c ~ 3) are 
preferred to explain the large values of ax observed at such scales. 



Although very recently deep Chandra of observations of Cen A (IHardcastle et al.ll2007l ) 
have provided us with spectra of some counterjet X-ray features, this is not the usual case 
for large scale counterjets. There are, however, some esti mates of the lower limi t of the 
X-ray brig htness asymmetry, like e. g., for the 3C 66B (> 25. Iliardcastle et al~ 2001 ). Pictor 
A (> 15, 



Wilson et all l200ll : 6 + ^ 2 Iliardcastle fc Crostonl hood ), or PKS 1127-145 (> 5, 



Siemiginowska et al.ll2002l ). Our model of AR or uniform radiating boundary layer tends to 
favor larger X-ray brightness asymmetries between the jet and the counterjet. The reason 
being that the X-ray emission of the jet is dominated by the synchrotron radiation from 
the monoenergetic, hard (7 = 10 8 ) electron component, while the counterjet emission is 
dominated by the EC radiation of the power-law part of the electron spectrum. 



5. Discussion 

SO02 chose the parameters of their radiating boundary layers such that they exhibit 
clearly the main characteristics of their model. We do not repeat here their discussion. 
Instead, we focus on the discussion of the choice of the new parameters of our model. 

SO02 anticipated that the form of the T(r) radial profile and of the related spatial 
variation of the acceleration efficiency can significantly influence the beaming pattern and 
the intensity of the boundary layer emission. We actually have confirmed this point by 
replacing the monotonically decaying profile of T(r) by an idealization of the Lorentz factor 
profile that results from AR boundaries. This results in a different relative normalization of 
the SSC and EC components presented in Figs. [5K 

The uniform layer model could be considered as a limiting case of AR layers when 
Tj } i — > T c . Thus, when looking at the plots of the jet-to-counterjet ratio (Fig. [8]) or the X-ray 
effective spectral index (Fig. [9]), any model including an AR layer with Iy G [10,15] (if 
T c = 10; black lines in Figs. [8] and [9]) or Ty e [3,4.5] (if r c = 3; grey lines in Figs. [8] and [9]), 
would display a graph located between the boundaries set by the solid and dashed lines of 
the corresponding figures. Since the area enclosed by these two line types (solid and dashed) 
is rather small, it turns out that our choice of the value of = 1.5r c does not significantly 
influence the results. More precisely, if for the considered values of T c (3 or 10) we would 
have picked up any other value of 6 [r c , 1.5r c ], neither the value of the jet-to-counterjet 
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brightness ratio, nor ax, nor even the spectrum would have changed appreciably. Since, 
keeping all other physical conditions fixed, a decrease of /3 h yields smaller values of Ty, the 
results we obtain shall be qualitatively valid both for magnetized or non-magnetized jets. 

The jet-to-counterjet effective X-ray spectral index asymmetry and its angular depen- 
dence (Fig. [9]) depends on which of the two spectral components at X-ray frequencies dom- 
inates: synchrotron radiation of the monoenergetic, high-energy electron component or the 
EC radiation from low energy electrons (7 ~ 100). In its turn, the dominance of any of these 
two processes depends on the exact shape of the electron spectrum at the highest energies, 
on the viewing angle, on the jet Lorentz factor and on T(r). Figure [9] serves as an example 
of the comparative behavior of different boundary layer models (in which all the parameters 
are the same except T(r)). 



5.1. Summary 



In this work, we compare the observational signatures imprinted by different kinematic 
models of shear layers of kiloparsec scale jets on the radiation of ultrarelativistic electrons 
accelerated at such boundary layers. The (simple) radiative model of a jet boundary layer 
follows that of SO02. Alternative (more sophisti cated) models to compute the b oundary 



layer radiation from multizone models (see, e.g., iMimica et all 120041 . 120051 . 120071 ) will be 
considered elsewhere. The most important difference among distinct kinematic models is 
the radial profile of the Lorentz factor T(r) across the shear layer. We have considered 
three different profiles. The first one is the usually assumed for large scale jets, namely, a 
monotonically d ecaying Lorentz f actor profi le (as in SO02). Th e second one is motivated 
by the results of lAloy fc Rezzollal (120061 ) and lMizuno et al.l (120081 ). namely, a non-monotonic 
profile where the Lorentz factor reaches a maximum at the limit between the jet and the 
external medium (AR layer). For completeness, the case of a uniform jet with a sharp edge 
is considered. The last case can be regarded as a limit of the AR profile when the maximum 
Lorentz factor in the boundary layer equals the Lorentz factor of the je t spi ne. Along the 



way, w e have conveniently developed the results of lAloy fc Rezzollal (120061 ) and lMizuno et al. 



(120081 ) to include the effects of dynamically important magnetic fields in the parameter range 
which is adequate for large scale jets. From such development we have obtained and idealized 
kinematic model of anomalous boundary layers in magnetized extragalactic jets. 

Not surprisingly, we find that the jet-to-counter jet brightness ratio at radio frequencies 
is larger for jets flanked by AR layers than for uniform or standard layers with the same jet 
core Lorentz factor. This results from the fact that the differences in the radial profile of 
the velocity field determine the beaming pattern of the boundary layer emission. Thereby, 
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this effect has to be taken into account when the jet bulk Lorentz factor is inferred from the 
jet-to-counterjet brightness asymmetry. 

The differences between AR and uniform boundary layer jets are small in the framework 
of our simple model. In practice, these small differences might render their observational 
distinction very difficult. Comparing these two models with standard radiating boundary 
layers the differences are much larger. Several independent clues can be used to distinguish 
(observationally) among them: 

1. If a sufficiently large number of TeV photons were detected, it would be possible 
to distinguish between jets with standard and anomalous layers by looking at the 
separation of the SSC and EC peaks. The separation between the two peaks is larger 
for jets with AR or uniform layers. 

2. A large jet-to-counter jet brightness ratio (about two orders of magnitude larger than 
in standard boundary layer jets) is expected for jets bounded by AR or uniform layers. 

3. For large viewing angles the effective X-ray spectral index is much larger for AR or 
uniform jets than for jets with standard boundary layers. 

4. There can exist an inversion of the jet and counterjet X-ray spectral indices for jets 
seen at large viewing angles. However, for most inclination angles and moderate bulk 
Lorentz factors of the jet core, the effective X-ray spectral index of the jet is much 
smaller than that of the counterjet. Indeed, if there is a hint on the bulk Lorentz 
factor of the jet (e.g., because superluminal proper motions are detected), the jet-to- 
counterjet spectral X-ray index ratio may tell us whether the radiating layer is uniform, 
standard or AR. We note that for large Lorentz factors of the jet core (r c ^ 8), the 
largest ratio a J x eff / eS corresponds to jets with uniform layers, whilst the smaller one 
should be identified with standard boundary layers. For more moderate core Lorentz 
factor, the ratio Q?xeft/ a x e ft * s l ar g er f° r jets flanked by AR layers than for jets with 
uniform or standard layers. In spite of these facts, the differences in the effective 
spectral X-ray index between models with uniform and with AR layers are rather 
small and, they might render their observational distinction very difficult unless a very 
careful analysis of the observational data is performed. 
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Fig. 1. — Left panel: Profiles of the Lorentz factor (solid line) and of the labora- 
tory frame magnetic field (dashed line) resulting from a prototypical Riemann problem 
yielding a ^SCS^ pattern. The initial left and right states are (p L , pi, v ™, T L , f3i) = 
(1(T 5 , 10" 7 , 0.9, 10, 0.1) and (p R , p R ,v^,T R , (3 R ) = (10" 6 , 10~ 2 , 0, 1, 0.1), respectively. Note 
that no boost is produced to the left of the contact discontinuity located at r ~ 1.147. 
Right panel: Same as the left panel, but for a prototype ^TZCS^ pattern. The initial left 
and right states are (p L , p L , v™, T L , /3 L ) = (10~ 5 , 10~ 7 , 0, 10, 100) and (p R , p R ,v^,T R , (3 R ) = 
(10~ 6 , 10~ 2 , 0, 1, 0.1), respectively. Note the boost (i.e., T* > TJ reached by the rarefied jet 
to the left of the contact discontinuity at r ~ 1.003. In both cases the initial discontinuity 
is set at r = 1 and the solution is shown after an arbitrary amount of time (we recall that 
the solution of the Riemann problem is self-similar). 
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Fig. 2. — Lorentz factor reached at the uniform state left to the contact discontinuity T* 
as a function of the normal velocity in the left state t>™. The right state is held fixed and 
is given by (p R , p R , v n , T R , j3 R ) = (1CT 6 , 10~ 2 ,0, 1,0.1). The left state has fixed values in the 
total pressure, the rest-mass density and the Lorentz Factor (p L ,p L ,T L ) = (10~ 5 , 10~ 7 , 10), 
while the magnetization fl h is varied for each curve as indicated by the different labels (see 
inset). The solid lines refer to the ^.TZCS^ pattern, and the dashed lines refer to a ^SCS^ 
pattern. The inset shows a zoom of the solutions for very small values of t>™, which are the 
most expected ones for standard jet flows. 
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Fig. 3. — Radial profile of the Lorentz factor across the jet showing the kinematic difference 
in the structure of an standard, monotonic shear layer (panel a) and an anomalous shear 
layer structure (panel b). 
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Fig. 4.- 



Schematic view of the jet model. 
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Fig. 5. — Predicted spectral energy distributions of the radiation generated in several bound- 
ary layers of a kpc-scale jet at 6* = 1° and 6* = 60° assuming 7 eq = 10 8 . The blue lines 
corresponds to the same monotonic shear layer profile as in SO02, the black lines to a uni- 
form boundary layer and the red lines to an anomalous one (see Sect. 13. 2p . The different 
contributions to the spectrum are labeled above the corresponding lines. The labels 1 and 2 
refer to the components associated to the power-law and to the monoenergetic parts of the 
electron distribution (Eq. [2D, respectively. 
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Fig. 6. — Predicted spectral energy distributions of the radiation generated by an AR bound- 
ary layer (T c = 10, Tj^ = 20) of a kpc-scale jet assuming 7 eq = 10 8 , for observing angles 
9* = 5°, 30°, 60° and 90°. Indices 1 and 2 denote the spectral components associated to 
the power-law isotropic energy distribution and to the monoenergetic part, respectively 
Absorption of VHE gamma-rays during the propagation to the observer is neglected. 
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Fig. 7. — Same as Figj6]but for a standard boundary layer (T c = 10, T^i = 1). 
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Fig. 8. — Jet-to-counterjet flux density ratio, S* (9*) / S* (n — 9*) at radio frequencies, as a 
function of the viewing angle 6*. Black and grey lines correspond to values T c = 10 and 
r c = 3, respectively -the values of r c are provided near the respective curves-. The black 
(grey) solid line corresponds to the radiation from an anomalous boundary shear layer with 
the radial profile of Eq. [14] and Tj i = 15 (Tj i = 4.5). The dashed-doted line corresponds to 
the radiation from a monotonic boundary layer with the radial profile of Eq. [14] and Tj^ = 1. 
The dashed line corresponds to the model with a uniform boundary layer with a Lorentz 
factor Tj 3 = 10. 
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Fig. 9. — Expected effective X-ray spectral index of a jet (j) and of a counterjet (cj) boundary 
layer emission, ax, e ff (#*), as a function of the viewing angle 9*. Solid, dash-dotted and dashed 
lines correspond to jets with anomalous, standard and uniform shear layers, respectively 
(Eq. HM . Black (grey) lines correspond to jets whose core moves with a Lorentz factor 
T c = 10 (r c = 3). The inset shows a zoom of the jet ax,efr for the smallest viewing angles 
9*. 



